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General remarks 

To achieve a quantitative understanding of the theory of 
strong interactions, QCD, at low energies, where the the- 
ory is non-perturbative, is one of the toughest challenges 
of todays particle physics. On the theoretical side in the 
last decades tremendous progress was achieved by two de- 
velopments, namely effective field theories (EFTs) as well 
as lattice gauge theory. In the former the symmetries of 
the fundamental theory are mapped on an effective the- 
ory formulated in terms of physical degrees of freedom 
(here it, K, r), D, D*, J/tp, ...) in a systematic and well- 
defined way based on a power counting in pertinent small 
parameter(s). In the latter the theory is solved on a dis- 
cretized space-time. Both approaches have a clear connec- 
tion to QCD, but they also have their drawbacks: Calcula- 
tions within lattice QCD are numerically very demanding 
and therefore a significant amount of expensive comput- 
ing time is necessary. For EFT calculations on the other 
hand need, when pushed to higher and higher accuracy, 
an increasing number of parameters — the so-called low 
energy constants (LECs) — needs to be determined, ei- 
ther by direct experimental input or using dispersion the- 
oryQ In recent years it became more and more clear that 
the future of the field lies in the combination of both these 
approaches, see e.g. (T). On the one hand, the LECs of 
the effective field theory related to the explicit symmetry 
breaking can be determined from varying the quark masses 
on the lattice - provided one is close enough to the chiral 
regime. In this way from a few lattice simulations in prin- 
ciple a large number of observables can be deduced, since 
the LECs typically contribute to various observables. Re- 
cent simulations allow to pin down the LECs related to the 
quark mass dependence of the pion mass and the pion de- 
cay constant to good precision, see e.g. (and references 
therein). Here, one can use the chiral perturbation theory 
expressions for these quantities to extract the LECs from 
the measured quark mass dependence. Another example of 
this interplay will be discussed below. On the other hand, in 
general lattice calculations are performed, at unphysically 
high quark masses, where computations are a lot cheaper. 
These can be connected to the physical parameter space us- 
ing extrapolation formulas derived using an appropriately 
tailored EFT. This strategy is illustrated in Fig. Q] EFTs 
also allow one to systematically explore the volume and 
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the lattice spacing dependences that arise naturally in any 
simulation. For a pedagogical introduction, see 0. 

In addition to the main lines from EFTs and lattice QCD, 
there is also shown the path of phenomenological mod- 
els. The large number of existing quark models fall into 
this class. The dashed arrows indicate that for those mod- 
els there is some connection to QCD — often the phrase 
"QCD inspired models" is used — and observables can be 
extracted, however, neither is the connection to QCD as 
rigorous as in the previous two examples nor is it normally 
possible to systematically improve the calculations. Still, 
phenomenological models like the non-relativistic quark 
model, can give important insights and are sometimes even 
used to estimate the free parameters of the EFTs. 

In addition to lattice QCD and EFTs, that allow one to 
calculate observables from first principles, in special cases 
also general theorems help to analyze data or to interpret 
results of calculations. One example of such a theorem 
applies, if the resonance of interest is located close to an 
s-wave threshold. Then the molecular admixture of this 
resonance by the corresponding continuum channel can 
be quantified model-independently. The resulting analysis 
method is outlined in the next section and applied to vari- 
ous states below. It should be clear that any field theoretical 
approach has to be consistent with such a general theorem. 
This holds especially for EFTs as well as lattice QCD. In 
this case one may use the mentioned theorem to interpret 
the results of some given calculation. An example of this is 
given below. 

In recent years a large number of narrow states was 
found in the charm-sector that do not fit into the conven- 
tional quark-antiquark picture. An equally large number of 
non-conventional explanations was proposed for those, in- 
cluding hybrids, glueballs, tetraquarks as well as hadronic 
molecules (forreviews, see e.g. Refs. !U|5]|6]|2][U), the last 
ones sometimes appearing as hadro-charmonia — bound 
systems of cc mesons and a light cloud (9). At present 
there are no model-independent methods available to dis- 
entangle the various scenarios mentioned, besides for the 
hadronic molecules as outlined in the previous paragraph. 
Thus, we will focus here on this kind of exotic states. 

Identifying molecular states 

Already 40 years ago Weinberg showed how the molec- 
ular content of the deuteron can be quantified fiOl . In 
Ref. [11] it was shown that this result can be generalized 
even to unstable states as long as their width is narrow and 
the inelastic thresholds are sufficiently far away. The cen- 
tral result of these studies is that the effective coupling con- 
stant of a resonance to a continuum channel in the s-wave 
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Figure 1 : Illustration of different methods to compare properties of QCD to experiment. 



can be written as 



mation on the true nature of the deuteron is hidden in the 
effective coupling constants in the sense discussed above. 



.9cff 



?f- = 4(m 1 +m 2 ) 2 A 2 v /2^ + , 

47T VP 

< 4(m 1+ m 2 ) 2 V2^+o(^p) (1) 

where e = m\ + m 2 — Mr denotes the binding energy 
of the resonance R, measured relative to the continuum 
threshold of particle 1 and 2, characterized here by their 
masses mi and m%. Further p, = m,\mil (mi + 7712) is the 
reduced mass and [3 denotes the range of forces. Clearly, 
the formulas given are useful only if the resonance pole 
is located very near a threshold, for only then the model- 
dependent corrections that involve the range of forces can 
safely be neglected. It should also be stressed that for all 
those states, where the nearby continuum channel is not 
in an s-wave, the scheme can not be applied. The most 
important parameter in Eq. (TTJ is A 2 . It gives the probabil- 
ity to find the molecular state in the physical state. Since 
the effective coupling is an observable, as will also be il- 
lustrated in the next paragraph, with Eq. ([1]) one is in the 
position to "measure" the amount of the molecular admix- 
ture in some physical state. Note, this implies especially 
that the structure information is hidden in the effective cou- 
pling constant, independent of the phenomenology used to 
introduce the pole(s). Therefore, a successful description 
of data using some model does not necessarily mean that 
it is based on correct assumptions. For example, if cer- 
tain spectra can be described using some isobar model — a 
model where all dynamics is parameterized by pole terms 
— it does not mean that all poles deduced are indeed gen- 
uine. The poles can as well parameterize singularities of 
the S'-matrix generated by hadron-hadron dynamics. A 
nice illustration of this fact is that one can very well study 
low energy nuclear physics in an effective theory with an 
explicit deuteron field 11121 . Also in this case, the infor- 



The 1^(4660) and 1^(4616) 

In order to illustrate how the method described in the 
previous paragraph works, in this section we will apply it 
to the F(4660). This state, discovered using ISR — fix- 
ing the quantum numbers to those of the photon, namely 
J PC = 1 — in the ip'irir final state |fl~3), is lo- 
cated very near the /o(980)^' threshold. It is therefore 
natural to ask, if it is an foip' molecule. An important 
feature of the F(4660) is that it was neither observed 
in e + e _ — + 7/si?. 7r+71 ' ~J/ip O, nor in the exclusive 
e+e~ -> DD, DD* , D* D* , DDn cross sections fl5), nor 



in the process e + e 



J/ijjD^D^ [16). In the f ip' 



molecular picture, these facts are easy to understand since 
the molecule would decay mainly through the decays of 
the unstable /o(980). However, they would severely chal- 
lenge other models which try to explain the Y (4660) as, for 
instance, a canonical 5 3 Si cc charmonium 1 17 1 or a baryo- 
nium 1 1 8 1 . 

As outlined above, for a given binding energy and given 
masses for the constituents, the effective coupling constant 
of a molecule can be calculated. Thus, it does no longer 
appear as a parameter in the fit. For the analysis we de- 
veloped an improved Flatte parameterization that allowed 
for a consistent inclusion of the spectral function of the /o 
via a dispersion integral — for details see Ref. [19]. As 
parameters in the fit to the data now only the mass of the Y 
and the overall normalization appeared. The fit gave 



M Y = (46651^) MeV 



(2) 



Clearly, by construction the mass of the Y (4660) is forced 
to stay below the f^ 1 threshold. The result of the fit, 
shown as the grey band, is compared to the data of Ref. lfl3l 
in Fig. 12 The fit is very good. Even the slight asymmetry 
visible in the data comes out naturally. To confirm that we 



were not misled by the fitting strategy we performed a sec- 
ond fit, where the effective coupling constant was allowed 
to float as well. Also this fit called for values of the cou- 
plings very near those extracted previously based on the 
molecular assumption, however, this more relaxed fit al- 
lowed for a larger range of masses for the Y. 
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Figure 2: Line shape for the Y (4660) under the assumption 
of a predominantly molecular nature. The best fit is shown 
as the solid line, the grey band reflects the 1-er uncertainty 
of the fit. The dashed line emerges as best fit, if also g c g is 
allowed to float. 

In Ref. COl it was argued that, if the F(4660) were 
an fotp' molecule, heavy quark spin symmetry predicts a 
pseudoscalar partner at My, = My — (M^i — M n < ) be- 
ing an /o?7c molecule with very similar properties. Espe- 
cially the line shape in the r]' c irir channel should be nearly 
identical to that shown in Fig. [2] The mass prediction, 
My = 4616j^g MeV, should hold with very high accuracy, 
since the spin symmetry breaking terms in case of hadronic 
molecules with a cc constituent are suppressed by an addi- 
tional factor I/ml — this emerges since those interactions 
need to be spin-dependent, connect colour neutral objects, 
and quark exchange is not possible. The width of the Y v 
into 7]' c irir is predicted to be T{r)' c Tnr) — (60 ± 30) MeV. 
An experimental confirmation of this prediction would be 
very desirable for it would provide further evidence for the 
picture sketched here. 

D* (2317) 

So far a method was described that allows one, under 
certain conditions, to spot molecules under the many states 
observed. However, the whole scheme was based on non- 
relativistic quantum mechanics. What we learn from the 
above is only additional information which states QCD 
forms, but no direct connection to QCD was employed. 
As outlined in the introduction, this is possible systemat- 
ically only using either lattice QCD or effective field theo- 
ries — more precisely chiral perturbation theory to control 
the dynamics of mesons formed of light quarks as well as 
heavy quark effective field theory to construct the interac- 
tions amongst heavy fields and, of course, suitable combi- 
nations thereof for heavy-light systems. 



In this section we will study resonances formed of a light 
and a heavy meson. However, no perturbative treatment is 
able to produce poles. Thus we need to leave the ground 
of rigorous effective field theories and add on top a resum- 
mation scheme. The power counting is done of the level 
of the effective potential so that the convergence of the se- 
ries has to be checked carefully for the resulting scattering 
amplitudes and bound state properties, much in the sense 
outlined by Lepage more than a decade ago EH . This 
resummation scheme at the same time has the beneficial 
side effect that all scattering amplitudes are unitary. Note, 
this procedure is completely analogous to what is routinely 
done in case of few-nucleon scattering — for a recent re- 
view see Ref. Il22ll . 

More precisely, we focus on the scattering of ir, K, and 
r) off D-mesons — see Refs. 1231 l24l . Our method is 
very similar to that of Ref. 11251 . however, for the first time 
higher order operators are studied completely and a com- 
parison to lattice data is performed. In addition the poles 
in the complex plane were identified revealing similar pat- 
terns in the heavy-light systems as observed previously for 
light resonances in Ref. |26|. 

It is especially important that a complete next-to- 
leading order calculation was performed which allowed for 
the first time for a reliable uncertainty estimate, which is, 
as we will see, especially important in case of the strong 
width of the _D* (2317), for here the quark model predic- 
tion is significantly lower OTl [32l than that derived within 
the molecular picture. The scalar charmed-strange me- 
son £>* (2317) is of particular interest since its observed 
mass is much lower than the expected value from the non- 
relativistic quark model [37|. Some ideas were proposed to 
shift the scalar charmed-strange meson mass down to the 
mass range of the _D* (2317). In ll38l . one loop corrections 
to the spin-dependent one-gluon exchange potential were 
considered. Another kind of modification is the mixing 
of the cs with the cqsq tetraquark [39|, or considering the 
coupling of the cs to hadronic channels, such as DK [ 40] . 
However, such kind of mass shift calculations are model- 
dependent as pointed out in Refs. 14T1I421 . A KD molecu- 
lar interpretation was proposed by Barnes et al. Il43l . which 
agrees with a lot of dynamical calculations [44, 27]. Other 
exotic explanations were also proposed, such as tetraquark 
state [45 1, and Dir atom [46|. Thus, in order to decide if 
an experiment can be decisive on the nature of this state, 
it is necessary that some measure is found to quantify the 
reliability of the results. 

Since the isoscalar Z)* (2317) is located below the DK 
threshold, the only open strong decay channel is D s tt°, 
which calls for isospin violation. In QCD there are two 
sources of isospin violation, namely quark mass differences 
and electromagnetism. Typically both contribute with sim- 
ilar strength, thus a scheme is necessary for the decays of 
the £>* (2317) that allows for a consistent inclusion of both 
type of effects. This framework is given by chiral pertur- 
bation theory. We performed a calculation up to next-to- 
leading order (NLO). It turned out, although essential for a 
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Figure 3: Goldstone boson scattering off D, D s -mesons in unitarized chiral perturbation theory. The dashed line shows 
the result of leading order chiral perturbation theory, while the grey band denotes the results of the full calculation. The 
lattice data are from Ref. 1291. 



quantitative understanding of the D°, D + mass difference, 
electromagnetic effects played a minor role in the strong 
decay of the D* (2317) — thus the decay was dominated 
by 7T° — r\ mixing, first shown in Ref. [27] to give a width 
of 8 keV, as well as the Z?-meson mass differences in the 
loops, studied previously in Refs. ||28l [25], which, when 
added in, give a width of 180 keV. On the contrary in case 
of a pure quark state the meson loop contribution is absent 
and thus a width of about 8 keV is predicted [32]. Similar 
values are expected for all other non-molecular scenarios. 
Thus, there is a direct experimental test available for the 
conjecture that the D* (2317) is a hadronic molecule. 

Unfortunately, the uncertainty of 110 keV turns out to 
be quite large. However, this uncertainty originates solely 
form the poorly determined isospin conserving interactions 
— to the order we are working all isospin breaking terms 
could be fixed from other sources. However, there are other 
tools to reduce the uncertainty on the isospin conserving 
part, namely by comparing to lattice QCD. For this com- 
parison scattering lengths (i.e the scattering amplitudes at 
threshold) are ideal: on the one hand they can be calculated 
on the lattice by systematic studies of the volume depen- 
dence of correlators [ 30 1 and on the other hand they can 
be straightforwardly calculated on the basis of (unitarized) 
chiral perturbation theory. The predictions of our calcula- 
tion are shown in Fig.[3]as a function of the pion mass. The 
dashed line shows the result of leading order chiral pertur- 



bation theory, while the grey band denotes the results of the 
full calculation. The lattice data are from Ref. [129]. Note 
that in the lattice calculations only the light quark mass was 
varied while the strange quark mass was kept fixed close to 
its physical value — this is why for vanishing m, only the 
scattering length for irD scattering vanishes, as demanded 
by the Goldstone theorem. As can be seen, in all chan- 
nels our calculation agrees very well with the lattice re- 
sults. The expansion appears to converge very well for irD 
scattering, while in case of DK and D S K scattering in the 
1=1 and 1 = 1/2 channels, respectively, up to 50% cor- 
rections emerged from higher orders. Although quite siz- 
able, even this change is within the expectations from the 
power counting as a result of the relatively large kaon mass. 
However, in case of KD scattering in the isoscalar chan- 
nel the change is truely dramatic — the scattering length 
even changes its sign. This change in sign, clearly a non- 
perturbative effect, is a consequence of the appearance of 
the D* (2317) as a KD bound state. In this context it is 
important to remember the theorem stated above that for 
bound states the effective coupling constant can be calcu- 
lated directly from masses, the binding energy and A 2 , the 
probability to find the molecular state in the physical state 
— c.f. Eq. ([T]i. This property translates directly into a pre- 
diction for the scattering length, if the binding energy of 
the bound state is small. In this case the scattering length 



may be written as 

^ -2A 2 /(1 + A 2 ) 

a ~ . (3) 

y/2mKe 

And indeed, the value we predict for the KD scattering 
length exactly agrees with the prediction from the gen- 
eral theorem in case of a purely molecular nature of the 
D* (2317) (A 2 = 0). This opens the opportunity to extract 
the nature of the Z?* (2317) directly from a lattice study: if 
this study were to find a scattering length in this channel of 
about —1 fm, it would prove the nature of the D* (2317) as 
a KD molecule. A larger absolute value of the scattering 
length is not allowed, while any smaller value would imme- 
diately quantify a non-molecular admixture — see Eq. ||3}. 
This equation also embodies the concept of universality for 
large scattering length, for a lucid review see ||47ll . 

Above we argued that lattice data can be used to de- 
termine the LECs of chiral perturbation theory. How this 
could work in the future can also be read off Fig. [3] At 
present both the theoretical predictions as well as the lat- 
tice data show a significant scatter. In addition, so far the 
lattice studies were performed only for one lattice spac- 
ing. Further investigations are necessary to quantify the 
systematics of the analysis and also adding more points for 
pion masses below 300 MeV. However, we hope that soon 
the lattice data/analysis will improve and then the parame- 
ters of our theoretical study can be constrained by a fit to 
those data (this calls for a complete theoretical calculation 
to NNLO, however, which still needs to be done). Once 
this is done, one can expect that the theoretical uncertainty 
of the prediction for the width of the £> s (2317) is reduced 
significantly as well. 

Discussion 

If mesonic bound states exist at all, shouldn't one ex- 
pect a large number of those in the spectrum? First of all 
one should stress that hadron dynamics normally produces 
sufficient attraction only in very few channels. This can 
also be seen, e.g., in Fig. [3] in all channels displayed, be- 
sides the KD 1 = channel, where the D* (2317) can be 
found, the interaction is repulsive. In addition, there may 
already exist various molecules in both, the light sector 
(here, e.g., the /q(980) is a prominent candidate for a KK 
molecule [33|, not only because of its mass close to the 
KK threshold but also because of its properties lfTTl[34ll ) as 
well as in the charm sector. Here the most prominent can- 
didate not discussed in this contribution is the X(3872) as 
a prime candidate for a D*D + h.c. molecule (see Ref. f35ll 
and references therein) or virtual state, possibly with some 
cc admixture (see Refs. [36] and references therein). In 
this contribution we outlined under what circumstances it 
is possible to get model-independent statements about the 
molecular nature of states — this is especially possible for 
states located very close to an ,s-wave threshold. However, 
models predict many more states. E.g. a brother of the 
F(4660) discussed above as a candidate for a foip' bound 



state, could be the F(4260) in Ref. O claimed to the a 
resonance of / (980) and J/ip. In Refs. EUED similar 
conclusions are drawn based on QCD sum rules. 

On the long run we will be able to address these issues, 
namely through the mentioned interplay of general theo- 
rems, lattice QCD as well as effective field theories. As 
was shown in Ref. [26], there are cases where it is possi- 
ble to move a resonance, by changing a QCD parameter — 
here the quark mass — from a position deep in the com- 
plex plane to a kinematic situation where the bound state 
analysis sketched above can be applied. The precondition 
for this is clearly to have a formalism which allows for a 
controlled quark mass dependence^ The parametric de- 
pendence of resonances on QCD parameters can also be 
checked in a study of the corresponding form factors l52l . 
This investigation provides an additional check of the sys- 
tematics of the approaches. In this sense we can hope to 
get a deeper understanding of QCD in the non-perturbative 
regime, once improved experimental data are available to 
check and refine the theoretical approaches. 
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